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Abstract 

In this article we study in detail the Q-vector dependence of the collective 

dynamics in simulated deeply supercooled SPC/E water. The evolution of the 

system has been followed for 250 ns at low T, allowing a clear identification of 

a two step relaxation process. We present evidence in favor of the use of the 

mode coupling theory for supercooled liquid as framework for the description 

of the slow a-relaxation dynamics in SPC/E water, notwithstanding the fact 

that the cage formation in this system is controlled by the formation of an 

open network of hydrogen bonds as opposed to packing constraints, as in the 

case of simple liquids. 

PACS numbers: 61.20.Ja, 64.70.Pf 
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I. INTRODUCTION 



The slow dynamics (a-relaxation) in supercooled molecular liquids and the glass- 
transition are two related topics which have received particular attention from the scien- 
tific community in the last years 0,0]. Significant progress have been made by a sinergetic 
approach based on theoretical, experimental and more recently on numerical work. Nowa- 
days, numerical simulations with realistic potentials describing the evolution of molecular 
system composed by thousand atoms, for time interval longer than 100 ns, are becoming 
feasible, allowing a closer check of the theoretical predictions and bridging the gap between 
experiments and theory. Such long simulations, although are suffering on the deficiency of 
the potential used compared to the real systems which they aim to simulate, offer a ideal 
reference for comparing with the novel theoretical predictions concerning the slow dynamics 
above the glass transition. 

In the case where the studied system is water, the interest in interpreting the molecu- 
lar dynamics in term of glass-transition concepts goes beyond the comparison between the 
simulated system and the theoretical predictions. Indeed, if the experimentally observed 
non-Arrhenius increase of the transport coefficients on supercooling can be explained 
within the same framework of simple supercooled liquids , the presence of hidden ther- 
modynamics anomalies at ambient pressure does not have to be invoked [^|Tni- the 
first paper of this series |Tl| we have tested Mode Coupling Theory (MCT) [T^ predictions 
for the correlation functions of single particle dynamics in water with corresponding quanti- 
ties calculated from Molecular Dynamics (MD) simulations, carried out for sufficiently long 
time (20 to 50 ns ) as to allow the slow dynamics to be observed. We tried to assess to 
what extent the MCT, which has been shown to describe simple liquids |jl3|,|l4|, is applicable 
also to the description of the single particle dynamics of (simulated) supercooled water, an 
hydrogen bonded liquid with strong non-isotropic interactions among molecules. As a result 
we found that the center of mass tagged particle dynamics can be qualitatively interpreted 
in terms of MCT. This result stimulated us to make the comparison more stringent and to 
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extend it to the collective center of mass dynamics. With one further year of computer time 
on three workstations working full-time, the limit of 250 ns of continuous simulation time 
has been reached, allowing the calculation of the collective properties, for which the noise 
level can not be reduced averaging over different molecules in the simulated system. 

Theoretical work based on the MCT has mainly focused on simple atomic liquids and on 
molecular liquids with spherical-symmetric interactions. Only recently the theory has been 
extended to treat the case of molecular systems with non-isotropic interaction potentials 
opening the way for a full and detailed quantitative comparison between theory and simula- 



tion/experiments for molecular systems [jl5l. Unfortunately, preliminary results at the ideal 
MCT level are only available for dipolar hard spheres |lT6| , p!7| , i.e. molecules with a simpler 
geometry than the water molecule geometry. Such results |16[ strongly supports the validity 
of general predictions based on the MCT for simple liquids. 

In this article we study in detail the Q-vector dependence of the collective dynamics in 
SPC/E water and present evidence in favor of a MCT description of the slow a- relaxation 
dynamics. The comparison is performed on a two level scale: (i) At qualitative level, we com- 
pare the MCT prediction for simple liquids (exemplified by the hard-sphere liquid within the 
Percus-Yevick approximation) with the center of mass dynamics in the case of supercooled 
SPC/E water. Although the comparison is by default qualitative, we think that the analogy 
in the Q dependence of all relevant parameters is particularly significant, (ii) At quantitative 
level, we compare the SPC/E slow dynamics with the dynamics predicted by MCT in the 
late /^-region; moreover we evaluate the von Schweidler exponents b and 7 governing the 
relaxation process (Sec.|I|) and we verify that they are related as MCT predicts. 

In both cases the agreement is striking and strongly supports the validity of MCT as 
the correct framework for interpreting the slow collective dynamics in simulated supercooled 
water. 
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II. MCT 



Several review papers present in detail the theory of mode coupling for supercooled 
liquids, in particular in its ideal formulation. A review on the predictions of the theory can 



be found in Ref. [11|. In this section we briefly report the main results which are relevant 
for the reading of the article. 

MCT aims at describing the slow dynamics in supercooled glass forming simple liquids 
or in molecular liquids with spherical-symmetric interactions between the molecules. It 
provides a description of the time evolution of density and current correlation functions 
(correlators in the following) in the time region where structural relaxation becomes the 
process which entirely controls the dynamics, i.e. for times longer than the microscopic 
characteristic time. MCT describes the slowing down of the structural relaxation, which is 
typical of supercooled liquids, including memory effects in a memory function which depends 
only on statical quantities (number density and structure factor S{Q)). 

In the ideal formulation the loss of correlation is ascribed completely to interaction 
between fluctuations of density pairs; all other channels for the decay of correlation, as 
for example the momentum modes, are completely neglected. The ideal MCT predicts a 
sharp liquid to glass transition, located at a finite critical temperature T^, associated with a 
power-law divergence of the correlation time. Tc is a purely kinetic transition temperature 
which describes the freezing of the molecules inside the cages and it does not deal with any 
thermodynamical anomalous behavior. The ideal-MCT can be considered as a first order 



approximation in a more complex scheme the kinetic transition is an artifact of the 
approximations involved and it disappears in the extended version of MCT, where couplings 
with particle momenta are taken into account. The power-law increase of the characteristic 
times predicted by the ideal MCT for T > Tc crosses to a different (unknown) relation, to 
describe the complex activated dynamical processes below Tc. The concept of Tc retains 
the meaning of cross-over temperature between two different dynamical behaviors and the 
prediction of the ideal MCT can be used to interpret the slow dynamics above T^ when 
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hopping effects are negligible. 

The ideal MCT predicts the existence of a two step relaxation scenario for the slow 
dynamics in supercooled simple liquids at temperatures close to Tc. The decay of the corre- 
lators, i.e. the loss of memory of the initial configuration, occurs in a two-step process P,p!9 



first the normalized correlator approaches a plateau value (pEA which is temperature- 
independent (non-ergodicity or Edwards-Anderson parameter); then decays from (pEA 
to zero (a-relaxation region). The short times region is strongly affected by the microscopic 
dynamics, which instead completely disappears in the long time region (a-relaxation). In 
the region close to the plateau (/3-relaxation region) (j){t) can be asymptotically expanded 
near the value (pEA- 

(pit) = (Pea + h^g{t/U) (1) 

where ta is the time scale which characterize the /3-region and a is the separation from 
the critical point (at fixed number density a = [T — T^/T^). The function g{t/t„) (also 
called /^-correlator) is solution of the asymptotic equation 

- 1 + \SLT[g\im + IM^)? = (2) 

where t = t/tcr is the scaled time, LT means the Laplace transform and z = zta- is the 
Laplace conjugate of t; A is the exponent parameter ||18| . 



It's easy to show that Eq. has a scale invariant solution; in the early /3-region, i.e. for 
t « t(j, g{t) has, at the leading order, the fractal behavior: 

^(t)=t-" + 0(n (3) 

In the late /3-region (t » t„) the predicted scaling law is 

git) = + Oii-") (4) 

The two scaling exponent a and b are not independent; they are both determined by the 
exponent parameter A through the relation: 



r(l-2a) r(i + 26) ^ ^ 

where F is the Euler's gamma function. The range of vahdity of these expansions for g{t) 
is hmited at small t by the microscopic dynamics time and at large t by the a-relaxation 
time. 

In the a-relaxation region another time scale appears, r; in terms of the rescaled time 
i = t/r the behavior of the correlator near the plateau is described in the leading order by 
a power law dominated by the exponent b. Including the next to leading order corrections 
the departure from (pEA in the early a-relaxation region is given by 

= ^EA - hl)i' + h2)t'' + 0{t^') (6) 

The amplitudes and h(^2) strongly depend on the physical features of the studied liquid 
and they have been explicitly calculated for a few models pO| , |2T|| . The a-relaxation time 
scale is a temperature dependent parameter which scales as the inverse of difFusivity: 

r(r) ~ |T-T,|-^ (7) 

with 

7 = l/2a + l/26 (8) 

The three scaling exponents a, b and 7 are not universal since they depend on the 
static structure factor. In the a-relaxation region the correlators obey the so-called time- 
temperature superposition principle, i.e. it is possible to scale the same correlator evaluated 
at different T on a single master curve, i.e. 

m = <p{t/T{T)) (9) 

For times much longer than r the mode-coupling dynamical equations can not be solved 
either analytically nor asymptotically. From numerical solutions developed for simple po- 
tentials it has been shown that a Kohlrausch- Williams- Watts (KKW) stretched exponential 
form, i.e. 

6 



(f){t) = AKexp[-{—)^'<] (10) 

can be used to empirically fit the last part of the (j)(t) decay. 

If we specialize the previous equations to the case in which is the intermediate scattering 
function F{Q, t), the self correlation function FseifiQ, t) or the transverse current correlation 
Jt{Q,t), the non-ergodicity factor 4>ea, the critical amplitudes h, /i(2), and the fitting 
parameters tk, Pk becomes Q-dependent quantities. Therefore the comparison between 
theory and experiments can be extended from the time dependence of and the temperature 
dependence of the correlation times to the Q-dependence of the above mentioned parameter. 

The MCT set of coupled integro-differential equations can be numerically solved once 
the static structure factor S{Q) is known. The problem has been solved for many different 
potentials: hard spheres [^^, soft spheres [^], Lennard- Jones [|^, binary mixtures of 



soft spheres and of Lennard- Jones [Q. For all these potentials the Q-dependence of 



the quantities appearing in the previous equations as well as the fitting parameters {tk and 
Pk) in Eq- shows pronounced oscillations in phase with S{Q). 

Moreover it has been shown theoretically [|],|19|,|18| that the large Q limit of the fitting 
parameters (3k coincides with the scaling exponent b. 

III. SIMULATION DATA 

The MD data analyzed in this article are a series of 250 ns long trajectories for a system 
of N=216 water molecules interacting with the SPC/E potential with periodic boundary 
conditions. SPC/E, one of the most common water- water potential, is a pair- wise additive 
potential designed to mimic the interaction between rigid water molecules via electrostatic 
and Lennard- Jones interactions. The simulation technique has been discussed in detail in 



Ref. and it is not repeated here. The only difference between the data analyzed in this 
article and the data analyzed in Ref. is the length of the simulation which now cover 
a time interval more than 5 times longer. We study seven different temperatures, ranging 
from T = 207K up to 285K. 



IV. ANALYSIS 



The basic quantity in the study of the center of mass collective dynamics in a liquid is 
the coherent intermediate scattering function, defined as 

1 ^ - 

FiQ.t) = ^ E e-'«[^^'=W-^^Wl (11) 
^W) j,k=i 

where N is the number of molecules in the system, Q is the wavevector and is the 
position of the center of mass of molecule i. F{Q, t) is the autocorrelation function of 
the space- Fourier transform of the density, thus giving information on the decay of density 
fluctuations at wavevector Q. 

The T and Q dependence of F{Q, t) are shown in Fig.|^ and ^ We show the T dependence 
at two different Q vectors, namely (i) the Q = 18nm~^ vector, corresponding to the position 
of the first pre-peak in the structure factor, Qfsdp (the analogous of the so-called first sharp 
diffraction peak in silica); (ii) the Q = 30nm^^ vector corresponding to the most intense 
peak in the structure factor, Qsqmax- In agreement with the general MCT predictions, all 
correlators show a two step decay process. For times smaller than 4ps we observe an oscil- 
latory decay process, from one to the non-ergodicity parameter. This time window contains 
all information about the microscopic dynamics and about the collective sound propagation 
in the Q range compatible with the simulation box size, a topic which have been studied 
in detail in several previous papers, independently from the MCT modelization and which 
will not be discussed in this article. We refer the interest reader to Refs. |]2^, where the 
frequency range above lOcm^^ ~ 27r/4 ps~^ has been studied in detail. The second decay 
process, from the non-ergodicity parameter to zero, has a strong T and Q dependence. This 
monotonic decay (the a-relaxation) becomes slower and slower on decreasing the tempera- 
ture. At the lowest simulated T, the a-relaxation is clearly separated from the microscopic 
time window and appears in a time region which was never explored before due to the long 
cpu time required to perform equilibrated molecular dynamics simulations. 

To quantify the Q and T dependence of the a-relaxation we fit all correlators for times 
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longer than 7 ps (to avoid the interference of the oscillatory sound modes) with a stretched 
exponential function (Eq.|TO|). For T higher that 230 K, the a relaxation time is smaller than 
20 ps and the a relaxation process, being superimposed to the collective sound modes, can 
be fitted with different parameters. Instead, F{Q,t) of the four lowest studied temperatures 
can be fitted unambiguously. The Q and T dependence of A^, tk and Pk, at the lowest 
four studied temperatures, are shown in Fig.||, together with S{Q). We note the presence 
of oscillations in all fitting parameters highly correlated with the oscillations in S{Q). 

To stress the similarity of our center of mass results with the prediction of MCT for 
spherical molecules we show in Fig.| Ak,tk and Pk evaluated from the numerical solution 
of the ideal MCT coupled integro-differential equations for the Percus-Yevick solution of the 
hard-sphere potential The MCT equations for the a-relaxation region have been solved 



using the same technique as in Ref. pO[. We note that both in the theoretical hard-sphere 
result and in the case of SPC/E center of mass dynamics, tk and f3 oscillate in phase with 
the structure factor. MCT also predicts that the limit of Pk at large Q coincides with the 
value of the exponent b in Eq.|^. From Fig.^we note that at large Q, Pk tends to the value 
0.5, the same value estimated previously on the basis of the analysis of the tagged particle 
motion for SPC/E water. 

According to MCT, close to Tc, the T dependence of tk{Q) at any fixed Q should follow 
Eq.|^. For the case b = 0.5, MCT predicts for 7 a value of 2.7. To compare with our data, 
we show in Fig. ^ tk{Q) /tk{Q fsdp) , i-e. imposing the equivalence of the scaled time at 
Q = Qfsdp- In a large Q range, i.e. Q < 40nm~^ all curves are nicely superimposed. We 
note that the longest relaxation time coincides with the position of the first sharp diffraction 
peak, as opposed to Qsqmax- Thus, the medium range order characteristic of network 
forming liquids is the most stable structure in the system. The enhanced stability at 
Q = Qfsdp can be predicted independently from MCT by calculating S{Q)/Q'^, i.e. on the 
basis of the so-called de Gennes narrowing. The inset in Fig.^ shows that the T dependance 
of the scaling factor tk{Qfsdp) is compatible with a power-law with exponent 7 = 2.7 and 
Tc = 2{)2K. We stress that the same value for 7 was found for the T dependence of the self 
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diffusion constant, but with Tc = 199K. The 2% difference in Tc is within the numerical 
error. The time-temperature prediction of MCT, i.e. the fact that all correlators have the 
same T-dependence of their relaxation time is consistent with the equal value of 7 found 
for the self and the collective dynamics in SPC/E water. As we discussed before in Ref. 
|TT|], the value of 7 calculated for SPC/E is different from the one obtained by fitting the 
experimental data for viscosity, diffusivity or reorientational NMR times — which run from 
1.6 to 2.4 on increasing the pressure. An extension of the present study at different isobars 
would be very valuable for clarifying the ability of the SPC/E potential to describe the 
dynamics of real water. 

We now turn to the SPC/E dynamics in the region where the value of the correlators is 
not very different from the non ergodicity factor, and where according to MCT the evolution 
of the correlators is controlled by Eq.p. According to MCT, curves at different T but at the 



same Q values fall on the same master curve. Following the procedure suggested in Ref. |T3 
we show in Fig. ^F{Q,t) vs. t/Ta{T) where is defined by F{Q,Ta) = 1/e, i.e. it is the 
time at which the correlation function has decayed to the 1/e value. From Fig. ^ it is clearly 
seen that, at all Q, curves for all T tends to sit on the same master curve, even the high T 
systems for which the a-region is hard to detect. To make the comparison with MCT more 
stringent, we compare the master curve designed by the envelope of the different samples 
with the universal master curve for the /3-region, i.e. the solution of Eq. As discussed 
above, the time dependence of g{t) is controlled only by the value of b, via Eq.|[ Thus we 
find that the same b value of 0.5, which was estimated from the study of self-dynamics in 
SPC/E water is able to rationalize both the limit value at large Q of as well as the time 
behavior of all correlation function in the /3-region. Unfortunately, no information can be 
obtained from our data concerning the exponent a, due to the superposition of the intense 
oscillations related to sound modes to the critical decay. 

As predicted by the theory, the range of validity of g{t) is different for different Q vectors, 
and appears to be larger for Q close to Qfsdp and smaller for Qsqmin = S7nm~^. The 
Q dependence of the correction to the master curves for large times (see Eq. ^ are also 
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predicted by MCT and can be compared with the center of mass collective dynamics in 
SPC/E water. To this aim we fit F{Q, t) at different Q at the lowest simulated temperature 
with the expression in Eq. ^ imposing b = 0.5, for all Q values. All correlator are fitted 
with the same b value. The quality of the fit is shown in Fig.|^, confirming that at all Q a 
rather good representation of the decay of correlation in the early a-region can be obtained 
in terms of Eq.^ Fig. |^ shows also that the contribution of the next-to-leading order 



correction is not at all negligible at some Q vectors |l20| , |21| . The result of the fit in the entire 
Q range are shown in in Fig.|[ Fig.^ shows the same quantities calculated theoretically for 
the hard-sphere case (PY approximation). As in the previous qualitative comparison, both 
h(i){Q) and h(i){Q) have oscillations in phase with S{Q). Results in Fig.H clearly highlight 
the need of performing an analysis in term of leading and next to leading expansion to detect 
the correct exponent b. In any case, as already discussed in Ref. F{Q,t) for large Q 
values (large meaning larger than Qsqmax) is not a good candidate for the identification of 
a von-Schweidler law. 



V. CONCLUSIONS 

In this article we have presented evidence in favor of a MCT description of the slow 
collective dynamics in deeply supercooled simulated SPC/E water. The presented data con- 
cerning the collective longitudinal dynamics, together with the data in Ref. [|lT| concerning 
the self motion, offer a complete picture of the dynamics of SPC/E water under deep super- 
cooling conditions. The evolution of the system has been studied for more than 250 ns at low 
T, allowing for the first time a clear identification of the two-step processes present in the 
correlators decay. The second decay process, clearly identified with the a-relaxation decay, 
is characterized by a highly non exponential behavior and by an apparent divergence of the 
characteristic time at a temperature around 200 K, i.e. about 50 K below the temperature 
of maximum density for SPC/E. An apparent divergence of the transport coefficient in real 
water, again at about 50 K below the temperature of maximum density, is also observed 
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experimentally [^,0 • 



The agreement between the data and MCT is striking both at qualitative and at quanti- 
tative level. At a qualitative level we have found that for both self and collective properties. 

(i) All correlators decay with a two step process which spreads over several time decades. 

(ii) The a-relaxation decay has an initial power-law behavior on leaving the plateau, whose 
range of validity is Q-dependent; the exponent in Eq. ^is 6 = 0.5; unfortunately an estimate 
of the exponent a is not feasible, due to the overlap in time with the sound modes, (iii) 
The long times part of the decay can be well fitted by a stretched exponential function, (iv) 
The Q-dependence of the parameters in the stretched exponential {i.e. amplitude, time and 
exponent) oscillates in phase with the static structure factor, in close analogy with the MCT 
predictions for simple liquids, (v) The Q-dependence of the von Schweidler amplitude and 
of the next-to-leading corrections are similar to the one predicted by MCT for a hard sphere 
system in shape and order of magnitude. 

At quantitative level we have found that (i) The stretching exponent (3 at large Q tends 
to 0.5, the same value as h as predicted, (ii) All characteristic times, both self and collective, 
satisfy the scaling law (Eq.|^ in a large Q range; moreover, the value of the exponent 7 is 
consistent with the value expected from the knowledge of the h exponent (Eq. ^ and it is 
independent from the correlator type, (iii) The correlator time dependence for values close 
to the non-ergodicity parameter is well described by the /5-correlator (Eq.Q). 

The analysis presented here and in Ref. represents an important step towards the 



understanding of slow structural relaxation in complex glass forming liquids, independently 
from the ability of SPC/E of mimicking real water. From the theoretical point of view, 
SPC/E is a non trivial molecular system, with molecules interacting via a strong anisotropic 
pair-wise additive interaction potential. Structural arrest in this system is not driven by 
packing constraints but by the formation of a strong network of tetrahedrally coordinated 
water molecules. Molecules in the cage are hindered in their translational and rotational 
motion by the presence of strong hydrogen bonds (energetic cages). The fact that MCT, 
which has been developed for simple liquids, succeeds well in describing the low T center of 
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mass dynamics in such a complicated molecular system, strongly supports the existence of 
a universality in the self and collective behavior of liquids under deep supercooling. Such 
universality has begun to emerge in the first applications to non-spherical molecules of the 
recently developed extension of the ideal MCT approach to treat angular correlators ||31|| ; 



this opens the way for an even more quantitative description of the SPC/E dynamics in 
a MCT framework. The series of simulation here analyzed can become a clean reference 
system to check the novel theoretical developements. 



The analysis of SPC/E water presented here and in Ref. |TT|] is also important in respect 
to the thermodynamic behavior of simulated water. We refer the interested reader to the 
recent book by Debenedetti and to Ref. |[T0|j9|,p^-|3^ . In this context, if SPC/E could 



be considered a sufficiently approximated modelling of real water, the data presented in 
this article would suggest that (if hopping processes did not intervene) at low pressure 
water would undergo a kinetic glass transition about 50 degrees below the temperature of 
maximum density, suggesting an interpretation of the so-called Angell temperature 
as the critical temperature of MCT P5|,|TT|. In this regard, the apparent power-law increase 



of transport coefficient in liquid water on supercooling is traced to the formation of cages 
and to the associated slow dynamics resulting from the presence of long living energetic 
cages. In other words the divergence of transport coefficients does not need to rely on a 
thermodynamical instability, either connected to the re-entrance of the gas-liquid spinodal 
or to the presence of a critical point at high pressure and low temperature l3^j3^ If SPC/E 
can represent the thermodynamic behavior of water sufficiently, we would conclude that 
at low pressure there is a continuous path connecting the hquid state to the low density 
amorphous ice. 
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FIGURES 

FIG. 1. F{QFSDP,t) (A) and FsQMAX,t) (B) vs time. QT = 207. K,aT = 210. K, 
T = 213.2 K, AT = 225.0 K,<T = 238.2i^, V ^ = 258.K, >T = 285. K Solid lines are fits 
with the KKW law (Eq.|l^) for times longer than 7 ps. Spherical average over all Q with the same 
modulus has been performed for all F{Q,t). 

FIG. 2. F{Q,t) at T = 207K. Q vectors are measured in units of Qo = O.lllnm^^. Solid lines 
are calculated according to EqJTo| 

FIG. 3. Fitting parameters of F{Q, t) at T = 207, T = 210, T = 215 and T = 225K according 
to the stretched exponential function Eq.^. Symbols as in Fig. [l|. 

FIG. 4. Fitting parameters of F{Q, t) for hard-spheres (FY) from the numerical solution of the 



MOT equations, according to the stretched exponential function Eq.lO 



FIG. 5. tk{Q) /tk{Qfsdp) as a function of Q. Note that in a large Q range, all curves 
can be scaled on one unique curve. The continuous line is the Q dependence of collective decay 
times based on the de Gennes narrowing hypothesis. Symbols in the inset show the scaling factor 
{tk{Qfsdp)) as a function of T. The full curve is the power law ~ |T — Td^'^, to highlight that 
the T dependence of the scaling coefficient is compatible with MOT predictions, with the predicted 
7 value. = 2Q2K 

FIG. 6. Master curve for F{QFSDP,t) (A), F{QsQMAX,t) (B), F{QsQMiN,t) (C). Symbols 
as in Fig.|^. Full line is the beta-correlator (the solution of the MOT equation ^ in the /3-region 
{g{t))). The dashed line indicates the F{Q,t) value chosen for scaling the different T. The three 
selected Q have been chosen to highlight the point concerning the Q dependence of the validity 
of the leading expansion t'^. The validity of the /3-correlator g{t) is a priori limited to the region 
where F{Q, t) — Fea is small. For clarity reason we have plotted g{t) in a larger range. 

FIG. 7. Fit according to Eq. ^ of the T = 207 K correlators limited to the early a-relaxation 
region {7ps { t { 800ps ), with h = 0.5. The fitting parameters are shown in Fig.^ 
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FIG. 8. Q dependence of the fitting parameters /i(2) to F{Q,t) at T = 207K. (See Eq.^ 
). and h(^2) have been multiplied by and t"^^ respectively to take into account 
the difference between i and real time. The arbitrary Q independent value of r has been 
chosen to bring and /i(2) in the same scale as the data of Fig.||, to make the qualitative 
comparison between the two set of data easier. 

FIG. 9. Q dependence of Fea, ^(i) and /i(2) for Hard-Spheres (PY) at the critical packing 
fraction r] = 0.51582. The PY solution for S{Q) is shown for reference. R is the hard-sphere radius. 



18 




t(ps) 




t(ps) 



FIG. 1. F. Sciortino et al 
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FIG. 2. F. Sciortino et al 
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FIG. 5. F. Sciortino et al 
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FIG. 7. F. Sciortino et al 
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FIG. 8. F. Sciortino et al 
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FIG. 9. F. Sciortino et al 
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